Current state of the modelling of photospheric 
activity 



A. F. Lanza 



Abstract I briefly review the current state of the modelHng of photospheric activity 
based on the high-precision optical light curves obtained with MOST, CoRoT, and 
Kepler. These models can be used to search for active longitudes where activity 
is preferentially concentrated, estimate the amplitude of stellar differential rotation, 
(-H I and look for short-term activity cycles as, e.g., in the case of CoRoT-2. In the case of 

a late-type star accompanied by a transiting hot Jupiter, the small light modulations 
observed during transits when a dark spot is occulted by the disc of the planet are 
also briefly considered. They can be used to derive information on individual active 
regions as well as on stellar rotation and the spin-orbit alignment of the system. 



1 Introduction 

In the Sun, we can study the photosphere in detail and map its brightness inho- 
mogeneities associated with convective motions, oscillations, and magnetic fields. 
In particular, magnetic fields give rise to solar active regions characterized by dark 
spots, bright faculae and an enhanced network. 

The transit of an active region across the solar disc as our star rotates produces 
a characteristic modulation of the optical flux with a rise when the region appears 
close to the approaching limb, a dip when it is closest to the central meridian, and 
another rise when it reaches the receding limb. The flux rises are associated with 
faculae that are brighter close to the limb and do not show an appreciable con- 
trast when they are close to the disc centre, while the flux dip is produced by the 
dark sunspots whose photometric effect is maximum when their projected area is 
maximum, i.e., when the active region is closest to the centre of the disc (cf., e.g., 
[Frohlich & Lean (2004)| ). 
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The collective behaviour of solar active regions shows a remarkable regularity. 
The total sunspot area oscillates with a mean period of 1 1 yr that is the main 
feature of the solar activity cycle. The latitudes where sunspot groups form migrate 
along the cycle, with the first groups appearing at ^ 30° — 35° and then migrating 
towards the equator as the cycle progresses. At the maximum of activity, sunspots 
form at a latitude of ^ 15° — 20° to reach a minimum latitude of 5° — 10° at the end 
of the cycle, when the sunspot groups associated with the next cycle begin to appear 
at high latitudes. 

The variation of the sunspot area is accompanied by a variation of the eigen- 
frequencies of the solar p-mode oscillations and of the amplitude and full width 
at half maximum (hereafter FWHM) of their peaks in the power spectrum. Since 
we wish to compare the solar behaviour with that of distant stars, we focus on 
low-degree modes (£ = 0, 1,2) that can be detected in disc-integrated observations. 
The variation of their eigenfrequencies along the previous solar cycle 23 and the 
onset of the present solar cycle 24 has been reported by jSalabert et al. (2009 )| , 
while I Jimenez-Reyes et al. (2003)) have presented the variation of the amplitude 
and FWHM of the mode peaks in the power spectrum. Note that the onset of the 
new cycle 24 has been apparent in the £~Q and 2 modes, while it has not been de- 
tected in the £ ^ I modes. This is related to the spatial distribution of the magnetic 
perturbations in the Sun because the ^ = 0,2 modes are more sensitive to the high- 
latitude perturbations, while the i = I modes respond mainly to the low-latitude 
perturbations. Since the magnetic fields are localized at high latitudes at the begin- 
ning of a new cycle, their effect can be detected in the ^ = 0,2 modes even before 
the first sunspot groups appear (cf. JSalabert et al. (2009)) for details). 

Analogous phenomena are expected in late-type stars with a level of magnetic 
activity comparable with the Sun. They have outer convection zones that host hy- 
dromagnetic dynamos generally leading to cycles ]Baliunas et al. (1995)| and dis- 
play convection-driven p-mode oscillations. Indeed, |Garcfa et al. (2010)) have de- 
tected a dependence of the eigenfrequencies and the mode amplitudes on the level 
of magnetic activity in the F-type star HD 49933, recently observed by the CoRoT 
(Convection, Rotation and Transits) space telescope. 

Taking the Sun as a template, to perform a comparative study of the impact of 
activity on p-mode oscillations in distant stars, we need to derive; a) the variation 
of the total spotted area vs. time, to trace stellar activity cycles; b) the latitude of 
the active belts vs. the phase of the activity cycle, assuming that there is a regular 
migration as in the Sun (see |Moss et al. (201 1)) for a detailed discussion of the pos- 
sible behaviours). In the next sections, we shall introduce the methods to map stellar 
photospheres that allow us to achieve these goals and present some application to 
the case of CoRoT targets. The advantages of stars accompanied by transiting hot 
Jupiters will be also emphasized. 
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2 Methods to map the photospheres of late-type stars 

To map stellar photospheres, we can apply methods based on spectroscopic and/or 
photometric data. Doppler Imaging (hereafter DI) is the most powerful technique 
based on spectroscopy, while photometric spot modelling is based on fitting wide- 
band light curves. 

DI is based on the line profile distortions induced by surface brightness in- 
homogeneities that can be resolved when the stellar rotation is the main source 
of line broadening |Vogt & Penrod (I983)| . In practice, this implies vsin/ > 10 — 
15 km s^' in a main-sequence star. In G or K-type dwarf stars, this corresponds to a 
maximum rotation period of ^ 4 — 5 days that implies an optical variability induced 
by magnetic activity of at least 0.02 — 0.05 mag. This may hamper the detection 
of p-mode oscillations as recently found by [Chaplin et al. (201 1)| by analysing a 



sample of Kepler targets (see also |Mosser et al. (2009a)| for the case of the active 
CoRoT target HD 175726). Another limitation of DI is the lack of sensitivity to a 
spot distribution that is uniform on length scales comparable with or larger than its 
spatial resolution. This poses a severe limitation on the measurement of the total 
spotted area in stars that are slowly rotating, as most of the asteroseismic late-type 
targets. Nevertheless, DI is capable of retrieving the latitude of starspots which is not 
possible with spot modelling, except in very favourable circumstances, and can be 
applied to spectropolarimetric observations of magnetic sensitive lines to map the 
surface magnetic fields and look for magnetic cycles (Zeeman Doppler Imaging, 
see, e.g., pares et al. (2009)| ). 

Spot modelling techniques are based on the reproduction of the rotational modu- 
lation of the optical flux induced by starspots. If there is a single spot on a rotating 
star, the time of light minimum can be used to derive its longitude and the maximum 
deficit of the flux gives a measure of the projected area of the spot. If the star is not 
viewed equator-on, the duration of the spot transit across the stellar disc can be used 
to estimate its latitude. As a matter of fact, several spots are usually present on the 
surface of an active star which makes the light curve inversion non-unique and often 
unstable because small changes in the data imply large changes in the derived spot 
positions and unprojected areas. In principle, a light curve is a one-dimensional data 
set, thus the derivation of a two-dimensional map of the photosphere without some 
additional assumptions is impossible. Nevertheless, the longitude of the spots and 
their total projected area are robust quantities that can usually be retrieved unam- 
biguously from a high-precision light curve. This has been demonstrated in the case 
of the Sun, for which we have spatially resolved maps of the photosphere, by e.g., 
[Lanza et al. (2007)[ . 

To alleviate the non-uniqueness and instability problems, two approaches have 
been applied. The first adopts a limited number of spots, usually up to 3 or 4, and 
fits a Hmited number of free parameters including the spot coordinates and areas 
(few-spot modelling). Several works have discussed the technique beginning from 
the earliest in the '70 and '80 [Torres & Ferraz Mello (l"973)l|Rodon6 et al. (1986)1 , 
to the first application to time series of data covering several consecutive rotations 
I Strassmeier & Bopp (1992)) , and finally to the time-series photometry obtained 
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by space-borne telescopes like MOST (the Microvariability and Oscillation of Star 
satellite, |Rucinski et al. (2004)[ ). 

The other approach to spot modelling consists in assuming a continuous dis- 
tribution of small spots, parameterized by a filling factor / giving the fraction of 
the surface occupied by the spots at each given location. To avoid the instability 
and non-uniqueness of the solution, the method needs a priori assumptions on the 
distribution of / to regularize the solution. The best approach, in the sense that it re- 
produces the distribution of sunspot groups in the Sun |Lanza et al. (2007)) , is based 
on the Maximum Entropy regularization (hereafter ME), where we look for the map 
with the maximum Shannon entropy among the infinite distributions of the filling 
factor that fit the light curve down to the precision of the measurements. From ME 
maps, the distribution of the filling factor vs. longitude and the variation of the total 
spotted area vs. time are derived. The information on spot latitude is usually very 
limited or absent in a light curve (cf. e.g., |Lanza et al. (2009a)) ), so we should not 
attempt to retrieve it from ME maps. 

Finally, we should mention the possibility of using the occultations of the 
starspots by a transiting planet as a high-pass spatial filter to map the fine structure 
of the spot distribution within the occulted band. The method has been successfully 
applied to, e.g., CoRoT-2 by |Silva-VaUo et al. (2010)||Silva-Valio & Lanza (20TT)) 
and provides information on the characteristics of individual spots as well as on 
stellar differential rotation and the spin-orbit alignment of the planetary system 
I INutzman et al. (2011)] . 



3 Results 

Long-term photometry allows us to derive stellar butterfly diagrams by tracing 
the systematic variation of the stellar rotation period as derived from the optical 
flux modulation vs. the time. The variation of the total spotted area can be mea- 
sured simultaneously by, e.g., the variation of the mean light level from one sea- 
son to the next, allowing us to trace the phase of the cycle and, in combination 
with the variation of the rotation period, the direction of migration of the spots 
[Messina & Guinan (2003)} . It is interesting to note that in some stars the activity 
belts may migrate polewards (see, e.g., ]Moss et al. (201 ij) ). In stars that rotate as 
slowly as the Sun, it is difficult to measure the rotation period using photospheric 
spots because they evolve on time scales of 10— 15 days in the Sun, i.e., signifi- 
cantly shorter than the rotation period itself. Chromospheric faculae live longer and 
are better tracers to look for the systematic variations associated with the solar cycle 



(see I Hempelmann & Donahue (1997) | and references therein). 

An estimate of the starspot temperatures based on multi-wavelength ground- 
based photometry has been presented by ]Amado & Zboril (2002)| . The high-precision 
space-borne photometry provided by MOST has allowed to measure the surface dif- 
ferential rotation of several late-type stars | Croll et al. (2006)] . Among the latest re- 



sults, those obtained by ]Mosser et al. (2009b) | are particularly interesting from the 
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point of view of asteroseismology because they refer to F-type main-sequence stars 
observed by CoRoT to look for p-mode oscillations. In addition to the spot lifetimes, 
they estimate the amplitude of the surface differential rotation and the inclination of 
the spin axis of the star to the line of sight with a typical uncertainty of 20° — 30°. 
Asteroseismology can do much better providing us with an inclination with an un- 
certainty of 5° — 10° in the case of such stars that rotate 4 — 5 times faster than the 
Sun |Gizon & Solanki (2003)) . In turn, this can be used to improve the estimate of 
the differential rotation from the spot modelling. 

The activity of several planet-hosting stars discovered by CoRoT has been inves- 
tigated. CoRoT-2 shows a short-term activity cycle with a period of about 29 days 
that could be analogous to the solar Rieger cycles [ [Oliver et al. (1998)l|Lou (2000)| 
[Zaqarashvili et al. (2010) or may be induced by some kind of star-planet interac- 



tion [Lanza et al. (2009a) 



[Lanza (201 1)| . Active longitudes have been detected and 



estimates of the surface differential rotation have been obtained for the other tar- 
gets CoRoT-4 ILanza et al. (2009b)| , CoRoT-7 ]Lanza etal. (2010)| , and CoRoT-6 
[ [Lanza etal. (201 1)| . 

The power of starspot occultation during transits in active stars accompanied by 
hot Jupiters has been shown by | Sanchis-Ojeda & Winn (201 1)| in the case of HAT- 
P-11 as observed by the Kepler space telescope. In addition to a quite accurate esti- 
mate of the inclination of the stellar spin, they have obtained a measurement of the 
true misalignment of the system by combining spot modelling with spectroscopic 
measurements of the Rossiter-McLaughlin effect. The latitude of the starspots has 
been constrained thanks to the large misalignment of the system. This offers a new 
opportunity to build stellar butterfly diagrams by means of the long-term Kepler 
photometry. 



4 Discussion and conclusions 

The recent detection of a dependence of the p-mode properties on the activity level 
in HD 49933, similar to that observed in the Sun, has revived the interest in studying 
stellar activity in late-type stars with outer convection zones. Those stars display p- 
mode oscillations excited by convection as well as activity phenomena produced by 
dynamo-generated magnetic fields. We have illustrated as the same high-precision 
photometric data collected to measure p-mode oscillations (and in several cases to 
search for planetary transits) can be used to model the distribution of starspots vs. 
longitude in a stellar photosphere and trace the variation of their area and migration 
along the activity cycle. This will allow us a better understanding of the impact of 
magnetic activity on p-mode oscillations in solar-like stars thanks to the long-term 
photometry provided by the CoRoT and Kepler missions (see the talk by S. Mathur). 
Spot modelling can significantly benefit from the determination of the stellar incli- 
nation made possible by asteroseismology. This in turn improves the estimation of 
the stellar differential rotation which is of fundamental importance to understand 
the rotational splitting of the p-mode frequencies and the operation of the stellar dy- 
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namo. Stars accompanied by a transiting planet offer new opportunities for the spot 
modelling thanks to the occultation of the spots by the planet and are also extremely 
interesting asteroseismic targets (see the talk by A. Moya). 
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